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, Abstract 

c ■ 

^ ■ Due to its relatively large eccentricity and proximity to the Sun, Mercury's orbital motion pro- 

\ vides one of the best solar-system tests of relativistic gravity. We emphasize the number of feasible 

ff^ ■ relativistic gravity tests that can be performed within the context of the parameterized weak field 

\ and slow motion approximation - a useful framework for testing modern gravitational theories in 

the solar system. We discuss a new approximation method, which includes two Eddington param- 
<^ \ eters (7,/?), proposed for construction of the relativistic equations of motion of extended bodies. 

QQ ' Within the present accuracy of radio measurements, we discuss the generalized Fermi-normal-like 

fS| \ proper reference frame which is defined in the immediate vicinity of the extended compact bodies. 

■ Based on the Hermean-centric equations of motion of the spacecraft around the planet Mercury, we 
\ suggest a new test of the Strong Equivalence Principle. The corresponding experiment could be per- 

■ formed with the future Mercury Orbiter mission scheduled by the European Space Agency (ESA) 
I for launch between 2006 and 2016. We discuss other relativistic effects including the perihelion 

■ advance, redshift and geodetic precession of the orbiter's orbital plane about Mercury. 
PACS number(s): 95.10Ce, 04.20.Jb, 04.20.Me, 04.70.-s 

U '. 

1 Introduction 

> ■ 

^ . It is well known that some gravitational effects depend not only on the first or second derivatives of the 

H ' gravitational potential, but also on the magnitude of this potential in the area where the gravitational 

■ - - ' experiments are performed. This is why, with the present level of experimental techniques, one should 

focus on locations with an intensive gravitational environment as the best place for conducting these 
studies. The planet Mercury plays a specific role in the history of modern gravitational physics. 
The successful explanation of the perihelion advance of its orbit has become one of the observational 
cornerstones for general relativity. Because of its proximity to the Sun, high eccentricity, and short 
orbital period, it offers a very interesting opportunity for the study of relativistic gravity. A spacecraft, 
placed in orbit about this planet, can provide the additional data necessary for dynamic tests of the 
principles of modern gravitational theories as well as the validation of the approximation methods 
used for analysis of the gravitational environment. 

In this paper we analyze the relativistic gravitational experiments for the Mercury Orbiter mission 
which has been included by the European Space Agency as a cornerstone mission under the Horizon 
2000 Plus program. The motivation for the research described here is to determine what scientific 
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information may be obtained during this mission, how accurate these measurements can be, and what 
will be the significance of the knowledge obtained. As it is known, there are three different types 
of measurements that are used in spacecraft navigation: radiometric (range and Doppler), very-long 
baseline interferometry (VLBI) and optical [Q]. In addition to navigation needs, the high precision 
Doppler, laser and radio range measurements of the velocity of and the distance to celestial bodies 
and spacecraft are presently the best ways to collect important information about relativistic grav- 
ity within the solar system. Combined with the technique of ground- and space-based VLBI, these 
methods provide us with a unique opportunity to explore the physical phenomena in our universe 
with very high precision. Most remarkable is that the accuracy of the modern VLBI observations 
is steadily increasing. Thus the delay residuals are presently of the order of 30-50 picoseconds (ps), 
which corresponds to an uncertainty in length of ~ 1 cm. Concerning the navigation of the interplan- 
etary spacecraft, the short arcs of spacecraft range and Doppler measurements, reduced with Earth 
orientation information referred to the International Earth Rotation Service' (lERS) celestial system, 
lead to a position determination in the extragalactic reference frame with an accuracy of order ~ 20 
milliarcsecond (mas). At the same time VLBI observations of the spacecraft with respect to extra- 
galactic radio-sources enable a direct measurement of one component of the spacecraft position in 
this extragalactic reference frame to an accuracy of about ~ 5 mas (||2|, [^). As a result, the use 
of such precise methods enables one to study the dynamics of celestial bodies and spacecraft with 
an unprecedented high accuracy. These data provide the necessary foundation for research of many 
scientific problems, such as: 

(i) . The construction of a dynamic inertial reference frame and a more precise definition of the orbital 

elements of the Sun, Earth, moon, planets and their satellites (Q, 

(ii) . The construction of a kinematic inertial reference frame, based on the observations of stars and 

quasars from spaceborne astronomical observatories (@, [^). 

(iii) . The construction of a precise ephemeris for the motion of bodies in the solar system to support 

reliable navigation in the solar system ( [|| , |^ ) . The construction of precise radio-star catalogs 
for the spacecraft astroorientation and navigation in outer space beyond the solar system. 

(iv) . The comparison of dynamic and kinematic inertial reference frames, based on the observations 

of spacecraft on the background of quasars, pulsars and radio-stars, as well as the verification of 
the zero-points of the coordinates in the inertial reference frame ([^], | p!o[ ], |pT|]). 

(v) . The experimental tests of the weak field and slow motion approximation (WFSMA) of modern 

theories of gravity ( Q , Q ) . 

These studies will enrich our knowledge about our universe, its cosmological evolution, and the 
behavior of stellar systems in general. By presenting this list of problems, we would like to emphasize 
the importance of the Mercury Orbiter mission not only for the gravitational theory, but for many 
of the fundamental problems stated above. It should be noted that, besides offering the opportunity 
to study the gravitational environment in Mercury's vicinity, this mission will also provide important 
data about the gravitational field of the Sun, Mercury's magnetosphere and its interaction with the 



solar plasma (see more detailed analysis of these problems in |14]), as well as enable one to verify the 
foundations of many recently proposed theories of relativistic reference frames ([|15|-[p|]). This mission 
provides a good opportunity for testing the principles of many modern approximation formalisms 
proposed for describing gravitational wave generation mechanisms (|jl2[, ||2^). It should be noted that 
a physically justified definition of the proper reference frame of the extended bodies {i.e. a complete 
multipolar solution of the gravitational N-body problem), in the considerably weak solar gravitational 
field, will allow us to better understand the physical processes in a stronger gravitational field regime. 
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This paper discusses the principles of a new approximation method developed for an astronom- 
ical relativistic N-body problem and analyzes the space gravitational experiments proposed for the 
future Mercury Orhiter mission. The outline of the paper is as follows. The next Section contains a 
brief introduction to the modern theoretical methods used to describe the motion of a system of N 
weakly interacting self-gravitating extended bodies and to analyze gravitational experiments within 
the solar system. We discuss problems associated with the traditional barycentric approach of the 
PPN formalism. In Section III we discuss the principles of a new iterative method for constructing 
the solution to the gravitational N-body problem in the WFSMA. In particular we discuss the physi- 
cal and mathematical properties of the relativistic astronomical reference frames. In Section IV this 
method is applied to derive gravitational field solutions and coordinate transformations. Our deriva- 
tions include two Eddington parameters (7, /?) which allow us to develop a new parametric theory 
of astronomical reference frames. Within the present accuracy of radio measurements, we discuss 
the generalized Fermi-normal-like proper reference frame, which is defined in the immediate vicinity 
of extended bodies. Thus, we have obtained the interesting result that although some terms in the 
Hermean-centric equations of motion of the spacecraft around the planet Mercury are zero for the 
case of general relativity, they may produce an observable effect in scalar-tensor theories. This allows 
us to propose a new test of the Strong Equivalence Principle (SEP), which is analyzed in Section V 
of the paper. Also in the fifth Section we discuss a number of relativistic gravitational experiments 
possible with the Mercury Orhiter. We present there both a quantitative and a qualitative analy- 
sis of the measurable effects such as Mercury's perihelion advance, the precession phenomena of the 
Hermean orbital plane, and the redshift experiment. In Section VI we present the conclusions and 
recommendations for future gravitational experiments with the Mercury Orbiter mission. In order to 
make access to the basic results of this paper easier, we present some calculations in the Appendices. 



2 Parametrized Post-Newtonian Gravity. 

Metric theories of gravity have a special position among all the other possible theoretical models. The 
reason is that, independent of the many different principles at their foundations, the gravitational field 
in these theories affects matter directly through the metric tensor of the Riemann space-time 
which is determined from the field equations of a particular theory of gravity. As a result, in contrast 
to Newtonian gravity, this tensor contains the properties of a particular gravitational theory as well 
as carrying the information about the gravitational field of the bodies. This property of the metric 
tensor enables one to analyze the motion of matter in one or another metric theory of gravity based 
only on the basic principles of modern theoretical physics. 

Within the accuracy of modern experimental techniques, the WFSMA provides a useful starting 
point for testing the predictions of different metric theories of gravity in the solar system. Following 
Fock ([^], perfect fluid is used most frequently as the model of matter distribution when 

describing the gravitational behavior of celestial bodies in this approximation. The density of the 
corresponding energy-momentum tensor T™-" is as follows^: 

(1) 

where p is mass density of the ideal fluid in coordinates of the co-moving reference frame, = dz^ /ds 
are the components of invariant four-velocity of a fluid element, and p{p) is the isentropic pressure 

*In this paper the notations are the same as in . In particular, the first three small latin letters a, h, c number the 
bodies and run from 1 to A^; the small latin letters k,l,...p run from to 3 and greek letters a,/?, 7,... run from 1 to 
3; the comma denotes a standard partial derivative and semicolon denotes a covariant derivative; repeat indices imply 
an Einstein rule of summation; round brackets surrounding indices denote symmetrization and square brackets denote 
anti-symmetrization. 
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connected with p by an equation of state. The quantity pYi is the density of internal energy of an ideal 
fluid. The definition of 11 is given by the equation based on the first law of thermodynamics (|l3|, 



11, m], my 

where p = y/—gpv? is the conserved mass density. Given the energy-momentum tensor, one may 
proceed to find the solutions of the gravitational field equations for a particular relativistic theory of 
gravity. The solution for an astronomical N-body problem is the one of most practical interest. In the 
following Sections we will discuss the properties of the solution of an isolated one-body problem as 
well as the features of construction of the general solution for the N-body problem in both barycentric 
and planeto-centric reference frames. 

2.1 An Isolated One-Body Problem. 

The solution for the isolated one-body problem in the WFSMA may be obtained from the linearized 
gravitational field equations of a particular theory under study. A perturbative gravitational field 
/i^", in this case, is characterized by the deviation of the density of the general Riemmanian metric 

tensor \/—gg^^ from the background space-time rjmn, which is considered to be a zero-th ordeiQ 
approximation for a series of successive iterations: yj—gg"^^ — \/^-r]i]mn = ^(o"' equivalently 



9mn — 



+ (3a) 



In order to accumulate the features of many modern metric theories of gravity in one theoretical 
scheme, as well as to create a versatile mechanism to plan gravitational experiments and to analyze the 
data obtained, Nordtvedt and Will have proposed the parameterized post-Newtonian (PPN) formalism 
(||2^-[^T[). This formalism allows one to describe the motion of celestial bodies for a wide class of 
metric theories of gravity within a common framework. The gravitational field in the PPN formalism 
is presumed to be generated by some isolated distribution of matter which is taken to be an ideal fluid 
Eq.(l). This field is represented by a sum of gravitational potentials with arbitrary coefficients: the 
PPN parameters. The two-parameter form of this tensor in four dimensions may be written as follows: 

4? = -2U + 2(/3 - r)C/2 + 2^ + 2r($2 - + (1 - 2zv)x,oo +0(0"'), (36) 
hSU = (27 + 2 - - T)Va + + t)Wo^ + 0(c-5), (3c) 

C = 2^?a/3(7 - r)U - 2TUap + 0(c-^), {3d) 

where rjmn is the Minkowski metric^ and the generalized gravitational potentials are given in Appendix 
A. 

Besides the two Eddington parameters (7,/?), Eq.(3) contains two other parameters v and r. The 
parameter v reflects the specific choice of the gauge conditions. For the standard PPN gauge it is given 
as z/ = ^, but for harmonic gauge conditions one should choose = 0. The parameter r describes a 
possible pre-existing anisotropy of space-time and corresponds to different spatial coordinates, which 
may be chosen for modelling the experimental situation. For example, the case r = corresponds 
to harmonic coordinates, while r = 1 corresponds to the standard (Schwarzschild) coordinates. A 

''For most of the non-radiative problems in solar system dynamics, this tensor usually is taken to be a flat Minkowski 
metric (@, @). 

®The following metric convention (H ) is used throughout as is the geometrical units h = c = G = 1. 
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particular metric theory of gravity in this framework with a specific coordinate gauge {u, r) may then 
be characterized by means of two PPN parameters (7,/?), which are uniquely prescribed for each 
particular theory under study. In the standard PPN gauge (i.e. in the case when 1/ = i,r = 0) these 
parameters have clear physical meaning. The parameter 7 represents the measure of the curvature 
of the space-time created by the unit rest mass; the parameter f3 is the measure of the non-linearity 
of the law of superposition of the gravitational fields in the theory of gravity (or the measure of the 
metricity). Note that general relativity, when analyzed in standard PPN gauge, gives: 7 = /3 = 1, 
whereas, for the Brans-Dicke theory, one has /3 = 1,7 = ^q^, where uj is an unspecified dimensionless 
parameter of the theory. 



The properties of an isolated one-body solution are well-known. It has been shown ( |13|| , |32], [p^ ]) 
that for an isolated distribution of matter in WFSMA there exist a set of inertial reference frames and 
ten integrals of motion corresponding to ten conservation laws. Therefore, it is possible to consistently 
define the multipole moments characterizing the body under study. For practical purposes one chooses 
the inertial reference frame located in the center of mass of an isolated distribution of matter. By 
performing a power expansion of the potentials in terms of spherical harmonics, one may obtain the 
post-Newtonian set of 'canonical' parameters (such as unperturbed irreducible mass /|^^ and current 

sj^^Q^ multipole moments^), generated by the inertially moving extended self-gravitating body (A) 
under consideration: 



^a(O) 



STF 



where f^"^ are the components of the symmetric density of the energy-momentum tensor of matter 
and gravitational field taken jointly. As a result, the corresponding gravitational field may be 
uniquely represented in the external domain as a functional depending on the set of these moments. 
Schematically this may be expressed as: 

where the functional dependence, in general, includes a non-local time dependence on the 'past' 
history^ of the moments [22|. However, by assuming that the internal processes in the body are 



adiabatic, one may neglect this non-local evolution. As a result, an external observer may uniquely 
establish the gravitational field of this body through determination of these multipole moments, for 
example, by studying the geodesic motion of the test particles in orbit around this distribution of 
matter 

It has been shown (||2^, (^S]) that for an isolated distribution of matter in WFSMA it is 

possible to consistently define the lowest conserved multipole moments, such as the rest mass of the 
body m^, it's center of mass z"^^, momentum and angular momentum 5""^. Thus the definitions 
for the mass rua and coordinates of the center of mass of the body z^^ in any inertial reference frame 
are given by the formulae (for more detailed analysis see ([13|, [34| and references therein): 



m« = / d^z- P(.'P), z:^{t) = ^ / d^z':: z':P\z'P). (5a) 



^To enable one to deal conveniently with sequences of many spatial indices, we shall use an abbreviated notation for 
'multi-indices' where an upper-case letter in curly brackets denotes a multi-index, while the corresponding lower-case 
letter denotes its number of indices, for example: {P} := /^i/^2 ■ • -Mpi ^{J} •S'^im-' - Pj • "^^^ explicit expression for 
the symmetric and trace-free (STF) part of the tensor r{p} is given in ([p6|-p^]). 

* Gravitational radiation problems are not within the scope of the present paper, and hence the set of multipole 
moments Eq.(4) are used for both tensor and scalar-tensor theories. 
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where the energy density iP^[z'p) of the matter and the gravitational field is given by: 

P\z'n =Pa{l+n + ^U- ^v^vf^ + 0(c-^)) , (56) 

with Pa being the conserved mass density defined by pa = Pa\/—gu^ and the velocity of the intrinsic 
motion of matter. In particular, the center of mass z"^ moves in space with a constant velocity along a 
straight line: Za^{t) = p'^ -t + k^-, where the constants = dz^^^/dt and /c" are the body's momentum 
and center of inertia, respectively. One may choose from the set of inertial reference frames the 
barycentric inertial one. In this frame the functions z^^ must equal zero for any moment of time. This 



condition may be satisfied by applying to the metric Eq.(3) the post-Galilean transformations [39| 
where the constant velocity and displacement of origin should be selected in such a way that p° and 
equal zero (for details see [^], pO[|). 

2.2 An Astronomical A^-Body System. 

By putting some restrictions on the shape and internal structure of the bodies, one can generalize the 
results presented above to the case of an isolated astronomical N-body system. Indeed, the assumption 
that the bodies posses only the lowest multipole moments (such as mass ma, intrinsic spin moment 
5"^ and the quadrupole moment Z"^) considerably simplifies the problem. The general solution with 
such an assumption is well-known (|^3|, see |^l| and references therein). The main properties of the 
solution Eq.(3) applied to such a case are well established and widely in use in modern ephemerides 
astronomy. In order to analyze the motion of bodies in the solar system barycentric reference frame, 
one may obtain the restricted Lagrangian function Lj^ describing the motion of N self-gravitating 
bodies ([|l3|, |^2[). Well within the accuracy necessary for our future discussion of the gravitational 
experiments on the Mercury Orbiter mission, this simplified function may be presented in the following 
form: 

N . N N / ^ 



1 3 1 

-(7 - T + -)Va^Va'' + - T + -)Va^Vb'' " (7 + 'r)nabxnab^iVa^Vb^ + T{nab^Va^f + 



riabX 



rab 



1 1 ^ / ^ mb\ 

(7 + 7j)l'a^ - (7 + l)^feM + nabxnab^^ ] + (/? " ^ " T^) E "^4 E " 

^abj ^ a ^b^a'^'^b' 

N N N 1 ^ 



-^EE E ^ambmc i^(^bc + ^ca) - — — +E"^'^'^(c (6) 

a bj^ac^a,b ^' ab ' ab' ac „ 

where is the isolated rest mass of a body (a), the vector r" is the barycentric radius- vector of 
this body, the vector r"j, = — is the vector directed from body (a) to body (fe), and the vector 
""afe ~ ^ab/'''ab is the usual notation for the unit vector along this direction. It should be noted that 
Eq.(6) does not depend on the parameter u, which confirms that this parameter is a gauge parameter 
only. The tensor /^^^ is the quadrupole moment of body (a) defined as: 

The tensor Sj^'^ is the body's intrinsic spin moment which is given as: 

sr = ^ [ d'z'^M^'j'MC - v^zn (8) 

''la J a 
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where vl^ is the velocity of the intrinsic motion of matter in body (a). Finahy, the quantity Ea is the 
body's gravitational binding energy: 



j3 2,'^(l^^"v Pa{^a)Pa{^ a' ) 
9m II a a \ ly III, 



The corresponding equations of motion of the planet in the solar system barycentric reference frame 
are given in Appendix B. It should be noted here that in the present numerical algorithms for celestial 



mechanics problems |13|], ||43|]-|^]) the bodies in the solar system are assumed to possess the lowest 



post-Newtonian multipole moments only. The corresponding barycentric inertial reference frame has 



been adopted for the fundamental planetary and lunar ephemeris |46|). Moreover, the coordinate 
time of the solar barycentric (harmonic) reference frame is the TDB time scale, which has been adopted 
by Fukushima [11| and Standish 



However, it turns out that generalization of the results obtained for the one-body problem, to the 
problem of motion of an arbitrary N-body system is not quite straightforward. Thus, we point out 
that taking into account the presence of any non- vanishing internal multipole moments of an extended 
body, significantly changes its equations of motion due to coupling of the intrinsic multipole moments 
of the body to the surrounding gravitational field. For example, for a neutral monopole test particle, 
the external gravitational field completely defines the fiducial geodesic world line which this test body 
follows (||l^, |23[|). On the other hand, the equations of motion for spinning bodies differ from these 
by additional terms due to the coupling of the body's spin to the external gravitational field (||48t|- 
||50|). As a result, one needs to present a post-Newtonian definition for the proper intrinsic multipole 
moments of the bodies in order to describe their interaction with the surrounding gravitational field 
and to obtain the corrections to the laws of motion and to the precession of extended bodies in this 



system ([|5l[-|56]). The fully relativistic definition of these moments may be given in the proper quasi- 
inertial reference frame only. Such a definition replaces Eq. (4) which was given in the rest-frame of the 
one-body problem^. In presenting these transformations one should also take into account that, due 
to the non-linear character of the gravitational interaction, these moments are expected to interact 
with external gravity, changing the state of motion of the body itself. Fock |^] was the first to notice 
that in order to find the solution of the global problem (the motion of the N-body system as a whole), 
the solution for the local gravitational problem (in the body's vicinity) is required. In addition, one 
must establish their correspondence by presenting the coordinate transformation by which the physical 
characteristics of motion and rotation are transformed from the coordinates of one reference frame to 



another. Thus, one must find the solutions to the three following problems (|17], |34|): 

(1) . The global problem: 

(i) . We must construct the asymptotically inertial reference frame. 

(ii) . We must find the barycentric inertial reference frame for the system under study. This is 

primarily a problem of describing the global translational motion of the bodies constituting 
the N-extended-body system {i.e. finding the geodesic structure of the space-time occupied 
by the whole system). 

(2) . The local problem: 



^Note that, due to the breaking of the symmetry of the total Riemann space-time by reaUzing the 3 + 1 spht 
these moments wiU not form tensor quantities with respect to general four-dimensional coordinate transformations. 
Instead, these quantities will behave as tensors under the sub-group of this total group of motion only, namely: the 
three-dimensional rotation. This is similar to the situation in classical electrodynamics, where electric E and magnetic 
H fields are not true vectors, but rather components of the 4x4 tensor of the electro-magnetic field Fmn = (i? ® H) 
0. 
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(i) . We must establish the properties of the gravitational environment in the proximity of each 

body in the system (i.e. finding the geodesic structure of the local region of the space-time 
in the body's gravitational domain). 

(ii) . We must construct the local effective rest frame of each body. 

(iii) . We must study the internal motion of matter inside the bodies as well as establish their 

explicit multipolar structure and rotational motion. 

(3). The theory of the reference frames: 



(i). We must find a way to describe the mutual physical cross-interpretation of the results 
obtained for the above two problems {i.e. the mapping of the space-time). 

Because the solutions to the first two problems will not be complete without presenting the rules of 
transforming between the global and the local reference frames chosen for such an analysis, the theory 
of the astronomical reference frames become inseparable from the problem of determination of the 
motion of the celestial bodies. ^From the other side, if one attempts to describe the global dynamics 
of the system of N arbitrarily shaped extended bodies, one will discover that even in WFSMA this 
solution will not be possible without an appropriate description of the gravitational environment in 
the immediate vicinity of the bodies ([0], l40[] ). 



3 Gravitational N-body Problem in the WFSMA. 

In this Section we will discuss the principles of a new iterative method for generating the solutions to an 
arbitrary N-body gravitational problem in the WFSMA of metric theories of gravity. In these theories 
one may choose any coordinates to describe the gravitational environment around the body under 
question. However non-optimal choice of these coordinates may cause unreasonable complications in 
the physical interpretations of the data obtained (see the related discussion in [^, |5^). Recently 
several different attempts have been made to improve the present solution to the N-body problem in 



WFSMA (see, for example, [|T5[-ll3], IQ, llH-[13])- Note that although these methods represent a 
significant improvement in our understanding of the general problem, the present situation requires a 
more physically adequate approach to account for the difficulties in describing the relativistic motion 
of a system of extended bodies. 

In the present paper we use an approach which has being developed to deal with the complica- 



tions introduced by using a standard PPN formalism |21|. This new method is based on the Fock- 
Chandrasekhar approach for dealing with extended and arbitrarily shaped celestial bodies and is 
parameterized with the two Eddington parameters (7,/3). This formalism is based upon the construc- 
tion of proper reference frames in the vicinities of each body in the system. Such frames are defined in 
the gravitational domain occupied by a particular body (b). One may expect that, in the immediate 
vicinity of this body, its proper gravitational field will dominate, while the existence of external gravity 
will manifest itself in the form of the tidal interaction only. Therefore, in the case of the WFSMA and 
in extreme proximity to the body under study, this proper reference frame should resemble an inertial 
frame and the solution Eq.(3) for an isolated one-body problem hmn! should adequately represent the 
physical situation. However, if one decides to perform a physical experiment at some distance from 
the world-tube of the body, one should consider the existence of external gravity as well. This is true 
because external gravity plays a more significant role at large distances from the body. 
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3.1 Definition of a Physically Adequate Proper Reference Frame. 

In order to construct a general solution for the N-body problem in a metric theory of gravity, let 
us make a few assumptions. First of all, we shall assume that the solution of the gravitational field 
equations /imn for an isolated unperturbed distribution of matter is known and it is given by the 
equations (3). With this assumption, we may construct the total solution of the global problem Qmn 
in an arbitrary reference frame as a formal tensorial sum of the background space-time metric r/mn, 
the unperturbed solutions /imn plus the gravitational interaction term TLmn- Thus, in the coordinates 
x'P = {x^,x'^) of the barycentric inertial reference frame, one may search for the desired total solution 
in the following form ||6^ : 



rimn{xn + E + -Hmnixn, (10) 

b=l 



where the coordinate transformation functions = y^(x^) are yet to be determined. The interaction 
term Timn will be discussed below. 

If the bodies in the system are compact and well separated, then, we may take into account that 
the mutual gravitational interaction between the bodies affects their distribution of matter through 
the metric tensor only. Therefore, without any loss of accuracy, we obtain the total energy-momentum 
tensor of the matter distribution in the system in the following form: 

Tmn{xn=j:§^^T,^dyB{xni (H) 
b=l 

where T^i is the energy-momentum tensoij^ of a body (b) as seen by a co-moving observer. 

One may establish the properties of the solution Eq.(lO) with respect to an arbitrary coordinate 
transformation simply by applying the basic rules of tensorial coordinate transformations. In partic- 
ular, in the coordinates = {yajUa) of arbitrary proper reference frame this tensor will take the 
following form: 

5mn(y2) = Q-mQ-n 9kl{x' {y^)) = q^q^ mix%y^a)) + 

+ g §§I^T(y'M) + ^§ (12) 

The expression for TmniVa) could be obtained analogously from that given by Eq.(ll). To complete 
the formulation of the perturbative scheme we need to introduce the procedure for constructing the 
solutions for the various unknown functions entering Eqs.(10)-(12), including the four functions of the 
coordinate transformations y^ = yj^iy^) and the interaction term Timn- 

^"As a partial result of the representation Eq.(ll) one can see that the Newtonian mass density pb of a particular body 
(&) is defined in a sense of a three-dimensional Dirac delta-function. Thus in the body's proper compact-support volume 
one will have: pb — mb5{yb), so that 

d^y'aPbiy'a) = mbSab, 



where Sab is the three-dimensional Kronekker symbol {Sb = Sab', Sab = 1 for a — h and — 0, for a 7^ b). Then in any 
proper reference frame the total density p of the whole N-body system will be given by the expression p(j/S) = 'Ylib P^'iVa)- 
This representation allows one to distinguish between the microscopic and macroscopic (or integral) descriptions of the 
physical processes and, hence, provides correct relativistic treatment of the problem of motion of an astronomical N-body 
system. 
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For the four functions of the coordinate transformations = ^^(2/^) ^se the general post- 
Newtonian coordinate transformations |21] which connect the coordinates (x^) of the barycentric 



inertial reference frame to those (y^) of a proper quasi-inertial reference frame of an arbitrary body 
(a). The general form of these relations may be given as follows: 

xO = + c-2/C.(yO, yl) + c-^£a(y°, yl) + ©(c"^), (13a) 

= Va + VaM) + c-^QZivl yl) + 0{c-^), (136) 

where the barycentric radius-vector of the body (a) in the coordinates of the proper reference frame 
is decomposed into Newtonian and post-Newtonian parts and is given as follows: 

M)) = VaoiVa) + ^ / d'y'a't^'iy'anQMyD + 0{c-% (14a) 



10 ^ 



m„ 



with 

/ d'y':P'{y'i) + 0{c-'), (146) 

J a 

where t^^{y'^) is the component of the conserved density of the energy-momentum tensor of matter, 
inertia and gravitational field taken jointly. The notation (..)^ applied to any function /(y^) means 
an averaging over body (a)'s volume as follows: 

(fiyl)) = f{yl) = — I d\:P,\y'nf{yly':)- (i5) 

\ /a rUa Ja 

In order to complete the formulation of the perturbative scheme, we need to introduce the procedure 
for constructing the solutions for the interaction term Timn and for the post-Newtonian transformation 
functions ICa, Q" and Ca which are still unspecified. The way to construct the solution for the interac- 
tion term Timn is quite straightforward; it is sufficient to require that the metric tensor in the form of 
Eq.(lO) or Eq.(12) will be the explicit solution of the gravitational field equations in the corresponding 
reference frame. These solutions are assumed to satisfy the covariant harmonical de Bonder gauge, 
which, for an arbitrary reference frame, may be written as follows: 



^tariyl) =0, (16a) 



where Pj^ is the covariant derivative with respect to the metric 'rf^^{y^) of the inertial Riemann- 
flat (RnmiiVstiyb)) ~ 0) space-time in these coordinates|^. For most of the practically interesting 
problems in the WFSMA in quasi- Cartesian coordinates this metric may be represented as the sum of 
two tensors: the Minkowski metric r]mn and the field of inertia (pmn- 

rLniyl) = = + <ttn{yl)- (166) 

Note that the term (pmn appears to be parameterized by the coordinate transformation functions 
lCa-,Ca and defined in Eq.(13); thus we have 4>mn{y^) = 4>mn[K-a-, ^^a-, Qa]i a formulation which will 
be referred to as the KCQ parameterization in the WFSMA. 



^^In Cartesian coordinates of the inertial Galilean reference frame {x^) the flat metric rfmn can be chosen as Tjmi = 
diag(l, —1, —1, —1), so that the Christoffel symbols Fmn^ = all vanish and conditions Eq.(16a) take the more familiar 
form of the harmonic conditions 

an(V=5<7™) =0, 

which are equivalent to setting = r = in the Eq.(3). 
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The search for the general solution for TCmnix^) is performed in a barycentric inertial reference 
frame (x^), which is singled out by the Fock-Sommerfeld's boundary conditions imposed on hmn and 
dkhmn (given by Eq.(lO)): 

r 

t -\ — = const, (17a) 
c 

where = —rj^i, x'^x'^. These conditions define the asymptotically Minkowskian space-time in a weak 
sense, consistent with the absence of any flux of gravitational radiation falling on the system from the 



external universe |12]. Moreover, one assumes that there exists such a quantity h'^'}^ = const (for the 
solar system this constant is of order 10~^) for which the condition 

hmnixn < hZT, im 

should be satisfied for each point (x) inside the system. Note that any distribution of matter is 
considered isolated, if conditions Eq.(17) are fulfilled in any inertial reference frame. 

The advantage of using the imposed conditions Eq.(17) is that it gives us the opportunity to 
determine the interaction term Ti-mnix^) in a unique way. It should be stressed that the corresponding 
solution gmn{x^) in the barycentric inertial reference frame resembles the form of the solution for an 
isolated one-body problem Eq.(3). The only change that should be made is to take into account the 
number of the bodies in the system: p — > J2bPbj where pij is the compact-support mass density of 
a body (b) in the system. However, both the interaction term Tlmn and the total solution for the 
metric tensor Qmn in the coordinates (y^) appear to be 'parameterized' by the arbitrary functions 
^aj^a, Qa- This result reflects the covariance of the gravitational field equations as well as the well 
defined transformation properties of the gauge conditions Eq.(16) used to derive the total solution. 
This arbitrariness suggests that one could choose any form of these functions in order to describe 
the dynamics of the extended bodies in the system. However, as we noted earlier, a disadvantageous 
choice of the proper reference frame (or, equivalently, the functions ICa,Ca and Q") might cause an 
unreasonable complication in the future physical interpretations of the results obtained. 



3.2 Principles of Construction of the Proper Reference Frame. 

In this Section we will present a way to find the transformation functions necessary for constructing 
a proper reference frame with well defined properties. As one can see from the expressions (12), 
in the WFSMA the main contribution to the geometrical properties of the proper reference frame 
in the body's immediate vicinity comes from its own gravitational field /imrT- Then, based on the 
Principle of Equivalence, the external gravitational influence should vanish at least to first order in 
the spatial coordinates (|^^, [^])- The proper reference frame, constructed this way, should resemble 
the properties of a quasi-inertial (or Lorentzian) reference frame and, as such, will be well suited 
for discussing the physical experiments. Note that the tensors hmn and TLmn represent the real 
gravitational field which no coordinate transformation can eliminate everywhere in the system. In the 
case of a massive monopole body, one can eliminate the influence of the external field on the body's 
world-line only. However, for an arbitrarily shaped extended body, the coupling of the body's intrinsic 
multipole moments to the surrounding gravitational field changes the physical picture significantly. 
This means that the definition of the proper reference frame for the extended body must take into 
account this non-linear gravitational coupling. 

In order to suggest the procedure for the choice of the coordinate transformations to the proper 
reference frame with well established properties, let us discuss the general structure of the solution 
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9mniya) given by expression (12). Thus, in the expressions for one may easily separate the four 
physically different terms. These terms are: 

(i) . The Riemann-flat contribution of the field of inertia (p^n given by expression (166). 

(ii) . The contribution of the body's own gravitational field h^n- 

(iii) . The term due to the non-linear interaction of the proper gravitational field with an external field. 

This contribution is due to the Newtonian potential and the potential $2 in the expressions (3). 
These interaction terms show up as the coupling of the body's intrinsic multipole moments with 
the external field. 

(iv) . The term describing the field of the external sources of gravity. This term comes from the 

transformed solutions hmn and the interaction term Timn- 

The first contribution depends on the external field in the gravitational domain occupied by the 
body (a) and appears to be 'parametrized' by the transformation functions Eq.(13). Note that for 
any choice of these functions, the obtained metric g'j^n^ by the way it was constructed, satisfies the 
gravitational field equations of the specific metric theory of gravity under study. Furthermore, based on 
the properties of the proper reference frame discussed above, one may expect that the functions /Ca, Ca 
and should form a background Riemann-flat inertial space-time t/^„ in this reference frame which 
will be tangent to the total gravitational field in the vicinity of body (a)'s world-line ja- Moreover, 
the difference of these fields should vanish to first order with respect to the spatial coordinates {i.e. 
the 'external' dipole moment equals zero [^3|). These conditions, applied to moving test particles, are 
known as Fermi conditions ( fs^ , |Q , |^ ) . We have extended the applicability of these conditions to 
the case of a system composed of N arbitrarily shaped extended celestial bodies. 

In order to obtain the functions ICa,Ca and for the coordinate transformation Eq.(13) we will 
introduce an iterative procedure which will be based on a multipole power expansion with respect 
to the unperturbed spherical harmonics. To demonstrate the use of these conditions, let us denote 
^mniUa) ^ ^he local gravitational field, i.e. the field which is formed from contributions (ii) and 
(iii) above. The metric tensor in the local region in this case can be represented by the expression: 
gmn\ya) = ^mri + H^^{yP). Then the generalized Fermi conditions in the local region of body (a) 
(or in the immediate vicinity of its world-line 7a) may be imposed on this local metric tensor by the 
following equations: 

7^7ti 7a 

hm rL(yD = ri(^)(yD , (i86) 

7^7a 7a 

where 7 is the world-line of the point of interest and the quantities rmn'^^(ya) are the Christoffel 
symbols calculated with respect to the local gravitational field gmn\ya)- Application of these condi- 
tions will determine the functions ICa, Ca, Qa which are as yet unknown. Moreover, this procedure 
will enable us to derive the second-order ordinary differential equations for the functions yaoiVa) and 
QaiUa^ 0), or, in other words, to determine the equations of the perturbed motion of the center of the 
local field in the vicinity of body (a). 

The relations Eq.(18) summarize our expectations based on the Equivalence Principle about the 
local gravitational environment of the self-gravitating bodies. By making use of these equations, we will 
be able to separate the local gravitational field from the external field in the immediate vicinity of the 
bodies. However, these conditions only allow us to determine the transformation functions for the free- 
falling massive monopoles (i.e. only up to the second order with respect to the spatial coordinates). 
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The transformation functions Eq.(13) in this case will depend only on the unperturbed contributions 
of the external gravitational potentials Uf, and and their first derivatives taken on the world- line of 
body (a). The results obtained will not account for the contribution of the multipolar interaction of the 
proper gravity with the external field in the volume of the extended body. This accuracy is sufficient 
for taking into account the terms describing the interaction of the intrinsic quadrupole moments of 
the bodies with the surrounding gravitational field, but some more general condition, in addition to 
Eq.(18), must be applied in order to account for the higher multipole structure of the bodies. 

The conditions Eq.(18), however, enable one to obtain the complete solution for the Newtonian 
function ICa- Functions Ca and may be defined up to the second order with respect to the 
spatial point separation, namely Ca,Qa ~ 0{\y^\^), so the arbitrariness of higher orders {k > 3) 
in the spatial point separation will remain in the transformation. In order to get the corrections to 
these functions up to A;*'' order (k > 3) with respect to the powers of the spatial coordinate y^, one 
should use conditions which contain the spatial derivatives of the metric tensor to order (k — 1). The 
mathematical methods of modern theoretical physics generally consider local geometrical quantities 
only and involve second order differential equations. These equations alone may not be very helpful 
for constructing the remaining terms in functions Ca, Qa up to the order k > 2. However, following 
Synge [^], one may apply additional geometrical constructions, such as properties of the Riemann 
tensor and the Fermi- Walker transport law ([^, p^-[|67| ). Another possibility is to postulate the 
existence of so called 'external multipole moments' (|15(|, |p^, [^] , ||37[|). However, those moments 
are defined through vacuum solutions of the Hilbert-Einstein field equations of general relativity in 
an inertial reference frame, while the influence of external sources of gravity are ignored. The fact 
of defining the moments in this way is essentially equivalent to defining the structure of the proper 
reference frame for the body under question. 

The most natural approach to define the desirable properties of the proper quasi-inertial reference 
frames for the system of extended and deformable bodies is to study the motion of this system in an 
arbitrary /C£Q-parametrized frame. There exist two different ways to do this, namely: (i) to study 
the infinitesimal motion of each element of the body, or (ii) to study the motion of a whole body with 
respect to an accelerated frame attached, say, to the center of inertia of the local fields of matter, 
inertia, and gravity. In our method we will use the second way and will study the dynamics of the 
body in its own reference frame. Our analysis will be directed toward finding the functions ICa,Ca 
and with the condition that the Riemann-flat inertial space-time ?7mn(ya) corresponding to these 
functions will be tangent to the total Riemann metric gmniUa) of the entire system in the body's 
vicinity. Physically, one expects that this inertial space-time will produce a 'fictitious' (or inertial) 
force with the density ficcQ acting on the body in its proper reference frame. At the same time, the 
body is under influence of the overall real force due to the local fields of matter and gravity with the 
density /q. Thus, the condition for finding the transformation functions ICa,Ca and is conceptually 
simple: the difference between these two densities T = fo — ftccQ should vanish after integration (or 
averaging) over the body's compact volume: 

6F= [ d^y'^f = I d^y'Jfo - f^ca) = 0. (19) 

Note that the notion of 'the center of mass' in this case loses its practical value, and one should 
substitute instead 'the local center of inertia'. Thus, the force fx.CQ. should provide the overall static 
equilibrium for the body under consideration in the local center of inertia, which is defined for all 
three fields present in the immediate vicinity of the body, namely: matter, inertia and gravity. Let 
us mention here that in practice it is not possible to separate these two forces /o and fK-CQ from 
each other. Fortunately, it is possible to obtain at once the difference between them T . This will 
considerably simplify the further analysis pi]]. 
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In order to construct the necessary solution for the functions /C^, Ca and in a way that will be 
valid for a wide class of metric theories of gravity, one must first analyze the conservation laws in an 
arbitrary /C£Q-parameterized reference frame. This could be done based on the conservation law for 
the density of the total energy-momentum tensor T"^"- of the whole isolated N-body system: 

V^f-"(y^) = 0, (20) 

where VJ^ is the covariant derivative with respect to total Riemann metric gmniVa) ™ these coordinates. 
Then, by using a standard technique of integration with Killing vectors, one will have to integrate 
this equation over the compact volume of the body (a), and one can obtain the equations of motion 
of the extended body (ll^, |24|, |p6[). Then the necessary conditions, equivalent to those of Eq.(18), 



may be formulated as the requirement that the translational motion of the extended bodies vanish in 
their own reference frames. This corresponds to the following conditions applied to the dipole mass 
moment rha = li^^: 

(Prha drha ^ ^ /oi ^ 

—-Z2 = -To" = "1-0 = 0, (21a) 

where the quantity is calculated based on the total energy-momentum tensor matter, inertia 
and gravitational field taken jointly (similar to condition of Eq.(4)). These conditions may also be 
presented in a different form. Indeed, if we require that the total momentum Pa of the local fields 
of matter, inertia and gravity in the vicinity of the extended body vanish, we will have the following 
physically equivalent condition: 

These conditions finalize the formulation of the basic principles of construction of the relativistic 
theory of celestial reference frames in the WFSMA. In the next Section we will present the results 
obtained when these principles are applied. 

4 Properties of the Proper Reference Frame. 

In this Section we will present the basic results obtained for the relativistic coordinate transformations 
in the WFSMA. We will present the results for the transformation functions Eq.(13) as well as the 
solution for the metric tensor in the coordinates of the proper reference frame. 

4.1 Coordinate Transformations and Metric Tensor. 

By taking into account all the conditions presented in the previous Section, one can obtain the set of 
differential equations on the transformation functions /Cq, Q2 and Ca- The solutions to these equations 
can be given as follows: 

/Ca(y°,yD = pt'{j2{Ub)^ - \va,Xo) - ^ao. " Va + 0(^-^)2/°, (22a) 

Qaivlva) = -tE (y^yf • {dpU,)^ - \yapyt{d''u,)^ + yZ{Ub)y 

b^a 

+yap fdt'[U:Xl + (7 + 1) E [(^f"^/^) + {d^'^U.vP^) ■ 

b^a 
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+ Kivl) + E Qa{L}{yA) ■ yi"-^ + o{Kt+') + o{c-%i, (226) 



/>3 



+ 



2 "cio '^do 



+ (7 + 1) E [(^[^F, 

6^a 



+ 



(7 + 1) - 2(7 + 1) E {vQ, - <(ya)_ 



+ 



(22c) 



Z>3 



where the dot over the function w'^^ denotes the regular time derivative. Note that in the case of the 
free-falling structureless test particle with conserved mass density given by /5a = i^aSiVa), the 
functions Eq.(22) correspond to the coordinate transformations to the proper reference frame defined 
on the geodesic world line of this particle. The time-dependent functions Q2{l} -^^a{z,} in the 
expressions Eq.(22) are the contributions coming from the higher multipoles {I > 3) (both mass and 
current induced ones) of the external gravitational field generated by the bodies (6 ^ a) in the system. 
These functions enable one to take into account the geometrical features of the proper reference 
frame with respect to three-dimensional spatial rotation. The form of these functions may be chosen 
arbitrarily. This freedom enables one to choose any coordinate dependence for the terms with Z > 3 in 
order to describe the motion of the highest monopoles. Moreover, one may show that, even though the 
total solution to the metric tensor gmn{x^) in the baryccntric inertial reference frame resembles the 
form of the one-body solution Eq.(3), if one will expresses this solution through the proper multipole 
moments of the bodies, it will contain the contributions from the functions Qa{x,].(ya) ^^d La{L}{y^)- 
Within the accuracy necessary for future analysis, we present the equations for both time-dependent 
functions and w'^^^ written with respect to time of the proper reference frame. Thus the 
Newtonian acceleration of body (a) with respect to the barycentric reference frame may be described 
as follows 



ao(ya) = -^' 



E^f d'y': Msf-^W^ + 0{c-% (23a) 



Whereas the post-Newtonian part of the acceleration may be represented by the expression 



^ b+a 



d 



d 



^(1 + 2^)v^v^ + ^ + (2/3 + 27 - 1)U)) + 

Z p 'la 
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l>3 



k 
l>3 



l>3 



Pa^T.QaiLM)d'^yi'' 
l>3 / 



9ya 



+ 0(c~ 



(236) 



For practical purposes one may find the value of the surface integrals in the expression Eq.(236) by 
performing an iteration procedure. It can be shown that the lowest multipole moments of the bodies 
will not contribute to this surface integration. However, the general results will fully depend on the 
non-linear interaction of the intrinsic multipole moments with the external gravity in the local region 
at the vicinity of the body under consideration. This additional iterative option will make all the 
results obtained with the proposed formalism easy to use in practical applications. 

The expressions Eq.(23) are the two parts of the force necessary to keep the body (a) in it's orbit 
(world tube) in the N-body system. These expressions are written in terms of proper time, and if one 
performs the coordinate transformation from the coordinates (y^) to those of (x^) for all the functions 
and potentials entering both equations (23), and takes into account the lowest intrinsic multipole 
moments of the bodies only, one obtains the simplified equations of motion for the extended bodies 
Eqs. (51-57) written in the coordinates (x^) of the barycentric inertial reference frame. 

The transformations Eqs.(13), (22)-(23) are the generalization of the Poincare' group for the case 



of motion of extended self-gravitating bodies. Note, that Chandrasekhar and Contopulos [39| obtained 
the coordinate transformations with the limitation of not violating the form-invariancy of the metric 
tensor. This ensured that the equations of motion would preserve their form. However, in our case, the 
transformation of coordinates to the proper reference frame connected with body (a), were obtained 
with the generalized Fermi conditions Eq.(18) and Eq.(21), and hence they changed the metric tensor 
considerably. The most notable contribution to the equations of motion of the test particle orbiting 
this body comes from the gravitational field of body (a) itself. The influence of the external sources 
of gravity presents itself in the form of the tidal terms only. This can be seen from the metric tensor 



in the coordinates (y^), which may be obtained in the following form: 



9miy'a) = 1-2U + 2W + yiiy^ ■ 7 V^.p A«ao " (2? 



1) "do /X '^"0/3 + 



+ 



d 



d 



k 



dy'a^' '^'dy^ 



1 + 



+2E [doLa{LM) + VaopdoQ^^^^iy'j] ■ + 0(|y;;|'=+^) + ©(c'^). 



l>3 



gUy^a) = 2(7 + 1) VaeV' + \l {MeS[ 



(24a) 
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M \k+l 



+0{\y\\'^+')+0{c-% 



(246) 



+ 



d 



d 



+ E [^o.xQ^a{L}{yl)i^+ripxQl{L}{yl)^\ ■yi'^^ + Om^-^^)+0{c-% (24c) 

;>3 Oya "Ua 

where the total gravitational potential U at the vicinity of body (a) is composed of the local Newtonian 
potential generated by the body (a) itself and the tidal gravitational potential produced by external 
sources of gravity: 



b bj^a 



(25) 



with the Newtonian acceleration of body (a) given by Eq.(23a). This potential is the solution of the 
Poisson equation in coordinates (y^), namely: 



(26a) 



By satisfying the requirement that the body is at rest in its proper reference frame, we derive the 
following results: 



d'y': PaiyTWivT) = I d'y': PaiyTWaiyT) = ^Ea, 
dU 



(266) 



'a'"- J a dyJr 



The quantity V°'{yP) in the expressions Eq.(24) is the total vector-potential produced by all the 
bodies in the system represented in the coordinates (yg) of the reference frame: 



+ v;- + 



bjta 



This potential satisfies the equation 



^ya (jya Ij 



(27) 



(28) 
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The form of the solution for 1/" Eq.(27) is chosen in such a way that it satisfies the Newtonian 
equation of continuity written in the proper reference frame: dop{yP) + dn[p{yP)v^^{yP)] = 0{c~^), 
which provides the following Newtonian-like relation between these quantities: 

dU dV 



dy'a dy^a ■ 

Another quantity we have introduced in the formulae Eq.(24) is W{y'^). This is the post-Newtonian 
contribution to the component 300 of the effective metric tensor in the coordinates (y^) of the proper 
reference frame and it is given as: 



with the expressions for the functions Wa and given as follows: 



(29) 



+ 



$:(2/?[/„(yS)C/ft(y2)- (37 + 1-2/3) / ^!^p„(yO,y;^-)t/,(yO,y^-) 



+ 



k 



d 



+ Y.Qa{L}{yl) / dh': Paiyly'^)^^ 

;>3 J a Oy^ 



ya ya 



\y'a - y'a'l 



d 



+ 0{\y:\''+')+0{c-% 
1 92 



(30a) 



c^a 



-(37 + 1-2/3) / ^!^p,(yO,y/.)^t/ft,(y°,y-)+ 

Wa ya I {,/ 



_ >{L} 

ya ya 



\y'a - y'^ 



+ 



+o(|yal'+') + o(0- 



(306) 



The functions Wa and W5 fully represent the non-linearity of the total post-Newtonian gravitational 
field in the proper reference frame. These functions contain the contributions of two sorts: (i) the 
gravitational field produced by the external (6 / a) bodies in the system, and (ii) the field of inertia 
caused by the accelerated (the terms with a^^) and non-geodesic motion (due to the coupling of the 
proper multipole moments of the body (a) with the external gravitational field and the self-action 
contributions both given by the terms with Qa{L} the expressions Eq.(30)) of the proper reference 
frame. 



4.2 The Fermi Normal Coordinates and the Equations of the Spacecraft Motion. 

The expressions Eq.(24) are the general two-parametric solution for the field equations of general 
relativity and scalar-tensor theories of gravity in WFSMA. These expressions satisfy the generalized 
Fermi conditions Eq.(18) at the immediate vicinity of body (a). This solution reflects the geometrical 
features of the proper reference frame with respect to the special properties of the motion of the A;-th 
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multipoles of the extended bodies. In order to find the unknown functions Qa{L}(ya) ™^ -^a{L}(y2) 
up to the k-th {k > 3) order, one should use the conditions which wih contain the spatial derivatives 
from the metric tensor to the 
conditions W 



— 1) order. Thus, following Synge [33| in addition to the Fermi 
one may apply the Fermi- Walker transport law 



Another possible method is to 
1). However, 



use the 'external' multipole moments as defined for the gravitational wave theory (|36| 
one can show that in WFSMA the functions Qa{L}^yV) ^"^^ ^a{L}{v^) may be chosen in such a way that 
the metric tensor Eq.(24) in a proper reference frame will take the form corresponding to any of these 
multipolar expansions. One can show that this happens because of the 'external' multipole moments 
in the transformation functions, and that this simply corresponds to the choice of coordinates of the 
reference frame with specific dynamical properties. Thus, for example, one can construct the proper 
reference frame for the Fermi normal coordinates |64]. Note that these conditions require corrections 
up to the third order with respect to the spatial coordinates in the transformation functions and 
>Ca The necessary corrections to the functions Eq.(22) have the following form|^: 



Va Va Va 



(31a) 



p ^ y Va Va Va ' ( 



d 



2{j + l)(d'^,Vhp 



-IV, 



10 u 



(316) 



Given this form for the corrections, the metric tensor Eq.(24) will take the form corresponding to 
generalized Fermi normal coordinates chosen in the proper reference frame: 



gUyP) = l-2U{yl) + 2Wa{yPa) 



+ 



(32a) 



dOaiVa) = 2(7 + 1) VaeV^iV^) + § ( 7 iVat^aaou - ri^uaaoa) + 



+(7 + 1) E [v.x{dl,V,') ^ - v.x{dUV,'') ^] ■ y^aV^'a + 0{K?) + 0(c-^), 
gaM) = ^o.p{l + 2^Ua{yl))+\lY. [ric.p{d%U,)^ + r,^,idlpUb)- 



-rifsJOM -r^aJd'M .yX + 0(|yal') + 0(c-"), 



(326) 



(32c) 



with the corresponding equation for a^^^ given by Eq.(23a). Thus we have obtained the form of the 
metric tensor in Fermi normal coordinates and the coordinate transformations leading to this form. 
These transformations are defined up to the third order with respect to the spatial coordinates. This 

With 52 is being the Kronekker symbol and rj^'^rjx,, = 
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accuracy is sufficient to analyze the motion of a spacecraft in orbit around Mercury consistent with 
parameterized relativistic gravity. 

We now obtain the equations of the spacecraft motion in a Hermean-centric reference frame. To do 
this, we consider a Riemann space-time whose metric coincides with the metric of N moving extended 
bodies. We study the motion of a point body in the neighborhood of body (a). The expression for 
the acceleration of the point body a^Q^ can be obtained in two ways: either by using the equations 
of geodesies of Riemann space-time du^ jds + T^^vJ^u^ = or by computing the acceleration of the 
center of mass of the extended body and then letting all quantities characterizing its internal structure 
and proper gravitational field tend to zero. In either case one obtains the same result. While the full 
derivation of the equations of motion is given in the Appendix C, we present here the restricted 
version of the equations (C4) which is consistent with the expected accuracy for ESA's Mercury 
Orbiter mission. This limited accuracy permits us to completely neglect contributions proportional to 
the spatial coordinates yj^. The planeto-centric equations of satellite motion around Mercury can be 
represented by a series in l/|y;,aol follows 




H 3 — I e\ 

Vbao 



Vbao y ytao 

(2/3 + ^7)r?-V + (/?-^KnV]'| + 0{\y'^\) + Oic-% (33) 



where 7^""^ = 77"'^ + '^Nj^^^N^^^ is the polarizing operator, subscript (a) denotes the planet Mercury 
and the post-Newtonian acceleration (^aO^o) '^^^ ^° gravitational field of Mercury only. This 



term is not new [TC] and it is given by the expression (C5a). Note that many relativistic terms in the 
geodetic equations have canceled out and, as a result, the equations of motion of a spacecraft around 
Mercury takes a very simple form Eq.(33). At this point we have all the necessary equations in order 
to discuss gravitational experiments with the future Mercury Orbiter mission. 



5 Gravitational Experiment for Post 2000 Missions 

Mercury is the closest to the Sun of all the planets of the terrestrial group and because of its unique 
location and orbital parameters, it is well suited to relativistic gravitational experiments. The short 
period of its solar orbit allows experiments over several orbital revolutions and its high eccentricity 
and inclination allow various effects to be well separated. In this Section we will discuss the possible 
gravitational experiments for the Mercury Orbiter mission. Analysis performed in this Section is 
directed towards the future mission, so we will show which relativistic effects may be measured and 
how accurately. 

By means of a topographic Legendre expansion complete through the second degree and order, the 



systematic error in Mercury radar ranging has been reduced significantly |71]. However, a Mercury 
Orbiter is required before significant improvements in relativity tests become possible. Currently, the 
precession rate of Mercury's perihelion, in excess of the 530 arcsec per century {" /cy) from planetary 
perturbations, is 43.13 "/cy with a realistic standard error of 0.14 "/cy After taking into account 
a small excess precession from solar oblateness, the later authors find that this result is consistent 
with general relativity. Pitjeva has obtained a similar result but with a smaller estimated error 
of 0.052 "/cy. Similarly, attempts to detect a time variation in the gravitational constant G using 
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Mercury's orbital motion have been unsuccessful, again consistent with general relativity. The current 
result I7| is G/G = (4.7 ± 4.7) x lO'^^ yr"!. 



Metric tests utilizing a Mercury Orbiter have been studied both at JPL and at the Joint Institute for 
Laboratory Astrophysics (JILA) and the University of Colorado. The JPL studies, conducted in the 
1970's, assumed that orbiter tracking could provide daily measurements (normal points) between the 
Earth and Mercury centers of mass with a 10 m standard error. A covariance analysis was performed 
utilizing a 16-parameter model consisting of six orbital elements for Mercury and Earth respectively, 
the relativistic gravity parameters (3 and 7, the solar quadrupole moment J20, and the conversion 
factor AU between unit distance (Astronomical Unit) and the distance in meters between Earth and 
Mercury. It was assumed that no other systematic effects were present, and that the normal-point 
residuals after removal of the 16-parameter model would be white and Gaussian. The total data 
interval, assumed equal to two years, corresponded to 730 measurements. Under the assumed random 
distribution of data, the error on the mean Earth-Mercury distance was 10/\/730 = 37 cm. The 
JPL studies showed, based on a covariance analysis, that the primary metric relativity result from a 
Mercury Orbiter mission would be the determination of the parameter 7, which describes the amount 
of spatial curvature caused by solar gravitation. The standard error was 0.0006, about a factor of 
two improvement over the Viking Lander determination. This accuracy reflected the effect of spatial 
curvature on the propagation of the ranging signal and also its effect on Mercury's orbit, in particular 
the precession of the perihelion. The error in the metric parameter /3 and the error in the solar J20 
were competitive with current results, but not significantly better. 



Within the last five years, a more detailed covariance analysis by the JILA group |74] assumed 
6 cm ranging accuracy over a data interval of two years, but with only 40 independent measurements 
of range. Unmodeled systematic errors were accounted for with a modified worst-case error analysis. 
Even so, the JILA group concluded that a two-order of magnitude improvement was possible in 
the perihelion advance, the relativistic time delay, and a possible time variation in the gravitational 
constant G as measured in atomic units. However, the particular orbit proposed by ESA for its 
2000 Plus mission was not analyzed. It is almost certain that the potential of the ESA mission lies 
somewhere between the rather pessimistic JPL error analysis and the JILA analysis of an orbiter 
mission more nearly optimized for relativity testing. 

In order to study the relativistic effects in the motion of the Mercury Orbiter satellite, we separate 
these effects into the three following groups: 

(i) . The effects due to Mercury's motion with respect to the solar system barycentric reference frame. 

(ii) . Effects in the satellite's motion with respect to the Hermean-centric reference frame. 

(iii) . Effects due to the dragging of the inertial frames. 

The effects of the first group are standard and all of them may be obtained directly from the 
Lagrangian function Eq.(6) or from the equations of motion Eqs.(i?l-i?7). The effects of the second 
group can be discussed based on the equations (33). And finally, the effects of the last group can be 
discussed based on the coordinate transformation rules given by Eq.(22). In the last case, however, 
we employ a simplified version of these transformations, due to the limited expected accuracy (~ 1 
m) of the Mercury ranging data. Thus, in the future discussion we will use the following expression 
for the temporal components: 

Ay^y^a) = yl + f fE — (i + (4^^ + /^^^"^^^^ 
1 



. _ „ . 2 

dt'-Vao,y^] +0{c-^), (34a) 
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where I^*^ represents the intrinsic quadrupole moments Eq.(7) of the bodies. The corresponding 
expression for the spatial components of the coordinate transformation is given by: 



Vbao 



with the precession angular velocity tensor given as follows: 

1, m;, 



K\yl) = E [(7 + )^^<A - (7 + i)5^<-S+ 



where S^^ is the intrinsic spin moments Eq.(8) of the bodies. 



(346) 



(34c) 



5.1 Mercury's Perihelion Advance. 

Based on Mercury's barycentric equations of motion one can study the phenomenon of Mercury's 
perihelion advance. The secular trend in Mercury's perihelion^ depends on the linear combination of 
the PPN parameters 7 and fi and the solar quadrupole coefficient J20 (|]l^], |75], |]76| ): 

'-Rm\2 J20nM 



^ = (2 + 27-/3)- 



+ 



3/ii0x2 



1 



p2 \2 



(3 cos iM - 1), "/cy 



(35a) 



«Af(l-e^^) ' 4VaM^ V^-^m; 

where aM^nMi^M and ejvf are the mean distance, mean motion, inclination and eccentricity of Mer- 
cury's orbit. The parameters /i© and Rq are the solar gravitational constant and radius respectively. 
For Mercury's orbital parameters one obtains: 



TT = 42" 98 



^(2 + 2, 



13) + 0.296 • J20 X 10"^ 



7cy 



(356) 



Thus, the accuracy of the relativity tests on the Mercury Orhiter mission will depend on our knowl- 
edge of the solar gravity field. The major source of uncertainty in these measurements is the solar 
quadrupole moment J20. As evidenced by the oblateness of the photosphere ||7^ and perturbations in 
frequencies of solar oscillations, the internal structure of the Sun is slightly aspherical. The amount of 
this asphericity is uncertain. It has been suggested that it could be significantly larger than calculated 
for a simply rotating star, and that the internal rotation rate varies with the solar cycle |78]. Solar 
oscillation data suggest that most of the Sun rotates slightly slower than the surface except possi- 
bly for a more rapidly rotating core |79]. An independent measurement of J20 performed with the 
Mercury Orhiter would provide a valuable direct confirmation of the indirect helioseismology value 
(2 ± 0.2) X 10"^. Further more, there are suggestions of a rapidly rotating core, but helioseismology 



determinations are limited by uncertainties at depths below 0.4 solar radii |8C]. 

The Mercury Orhiter will help us understand this asphericity and independently will enable us to 
gain some important data on the properties of the solar interior and the features of it's rotational 
motion. Preliminary analysis of a Mercury Orhiter mission suggests that J20 would be measurable to 
at best ~ 10~^ [74| or about 1% of the expected J20 value. This should be compared with the present 
10% solar oscillation determination [f77| . 

It should be noted that the Mercury Orhiter itself, being placed in orbit around Mercury, will experience the 
phenomenon of periapse advance as well. However, we expect that uncertainties in Mercury's gravity field will mask the 
relativistic precession, at least at the level of interest for ruling out alternative gravitational theories. 
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5.2 The Redshift Experiment. 

Another important experiment that could be performed on a Mercury Orbiter mission is a test of 
the solar gravitational redshift. This would require a stable frequency standard to be flown on the 
spacecraft. The experiment would provide a fundamental test of the theory of general relativity and 
the Equivalence Principle upon which it and other metric theories of gravity are based pi]]). 
Because in general relativity the gravitational redshift of an oscillator or clock depends upon its 
proximity to a massive body (or more precisely the size of the Newtonian potential at its location), 
a frequency standard at the location of Mercury would experience a large, measurable redshift due 
to the Sun. With the result for the function /Ca given by Eq.(22a) and Eq.(34a) in hand, one can 
obtain the corresponding Newtonian proper frequency variation between the barycentric standard of 
time and that of the satellite (the terms with the magnitude up to 10~^^), given as: 

= 1 + + ^ + + mf - ^(RMVio)) + 0{c-% (36a) 

dy(o) Rm y{o) 2c^ R%j ^ ' 



where {y^Qy 2^(0)) four-coordinates of the spacecraft in the Hermean-centric reference frame and 

{/(o) is the spacecraft orbital velocity. One can see that the eccentricity of Mercury's orbit would be 
highly effective in varying the solar potential at the clock, thereby producing a distinguishing signature 
in the redshift. The anticipated frequency variation between perihelion and aphelion is to first-order 
in the eccentricity: 

f^) =^J^. (366) 

This contribution is quite considerable and is calculated to be ((J///o)eM = ^-^ ^ 10^^. Its magnitude, 
for instance, at a radio-wave length Aq = 3 cm (/o = 10 GHz) is {5f)eM — HO Hz. We would 
also benefit from the short orbital period of Mercury, which would permit the redshift signature of 
the Sun to be measured several times over the duration of the mission. If the spacecraft tracking 
and modelling are of sufficient precision to determine the spacecraft position relative to the sun to 
100m (a conservative estimate) then a frequency standard with 10~^^ fractional frequency stability 
5f / f = 10^^^ would be able to measure the redshift to 1 part in 10^ or better. This stability is within 
the capability of proposed spaceborne trapped-ion |82] or H-maser clocks ([^, |84]). 



5.3 The SEP violation effect. 



Besides the Nordtvedt effect (for more details see |14| and [^]), there exists an interesting possibility 



for testing the SEP violation effect by studying spacecraft motion in orbit around Mercury. The 
corresponding equation of motion is given by Eq.(33). As one can see, the two terms in the second 
line of this equation vanish for general relativity, but for scalar-tensor theories, they become responsible 
for small deviations of the spacecraft motion from the fiducial geodesic. Both of these effects, if they 
exist, are due to non- linear coupling of the gravitational field of Mercury to external gravity. They 
come from the expression for Wa given by Eq.(28o), which is the local post-Newtonian contribution 
to the qqq component of the metric tensor in the proper reference frame. 

The first of these terms may be interpreted as a dependence of the locally measured gravitational 
constant on the external gravitational environment and may be expressed in the vicinity of body (a) 
as follows: 

Ga = Gq ^ ^ „ , ^ 



l-(4/3-37-l)E— 1- (37) 
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In the case of a satellite around Mercury, the main contribution to this effect comes from the Sunf^. 
Because of the high eccentricity of Mercury's orbit, the periodic changing of the sun's /oca/ gravitational 
potential may produce an observable effect, which can be modeled by a periodic time variation in the 
effective local gravitational constant: 



^/ifl Q 1^^ /^Q iiceMsin0(t) 2 /goN 

(4/3-37-1) Ti 2^{l + CM cos (f>{t) ) , (38) 



aA/(l-e^j)J aj\/(l 



T J period 

which gives the following estimate for this effect on Mercury's orbit: 

- (4/?-37- 1) X 1.52 X 10"^sin(/)(t) yr'^ (39) 

- Gr J period 

Note that this effect Eq.(39) is fundamentally different from that introduced by Dirac's hypothesis 



of possible time dependence of the gravitational constant ||7^. As one can see from expression (39), 
the characteristic time in this case is Mercury's siderial period. This short period may be considered as 
an advantage from the experimental point of view. In addition, the results of the redshift experiment 
could help in confident studies of this effect. Recently a different combination of the post-Newtonian 
parameters in the Nordtvedt effect, rj = 4/3 — 7 — 3, was measured at 77 < 10~^ [^. This means 
that, in order to obtain the comparable accuracy for the combination of parameters Eq.(39), one 
should perform the Mercury gravimetric measurements on the level no less precise than [G/G]^^^^^^ ~ 
10^^^ yr^^- Recently a group at the University of Colorado has analyzed a number of gravitational 



experiments possible with future Mercury missions [74|. Using a modified worst case error analysis, 
this group suggests that after one year of ranging between Earth and Mercury (and assuming a 6 
cm rms error), the fractional accuracy of determination of the sun's gravitational constant ttiqG is 
expected to be of order ~ 2.1 x 10~^^. Moreover, even higher accuracy could be achieved with a 
Mercury lander as proposed by ||7^ . This suggests that the experiment for determination of the effect 
Eq.(39) may be feasible with the Mercury Or6iier mission. 

Another interesting effect on the satellite's orbit may be derived from the Eq.(33) in the form of 
following acceleration term: 

<5a(o)s£P = 2(/? - l)!^^^(iVM - n(niVM)), (40) 

where Rm is Mercury's heliocentric radius- vector and Nm is the unit vector along this direction. 
This effect is very small for the orbit proposed for ESA's Mercury Orbiter mission. However, one 
can show that there exist two resonant orbits for a satellite around Mercury, either with the orbital 
frequency a;(o) equal to Mercury's siderial frequency ljm- '^(0) ~ or at one third of this frequency 
u;(o) ~ cjaz/S. For these resonant orbits, the corresponding experiment could provide an independent 
direct test of the parameter /3. 



5.4 The Precession Phenomena. 

In addition to the perihelion advance, while constructing the Hermean proper reference frame, one 
should take into account several precession phenomena included in the transformation function 
and associated with the angular momentum of the bodies. As one may see directly from Eq.(226) and 

^''Note that this combination of PPN parameters differs from that for a similar effect presented in jl^]. The reason 
for this is that, in our case the transformations in the form of Eq.(22) enable us to define the physically adequate 
transformation rules for the metric tensor between the barycentric and planeto-centric reference frames and, hence, to 
obtain the correct and complete equations of geodesic motion Eq.(33) in the Hermean-centric reference frame. 
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Eq.(346), besides the obvious special relativistic contributions, the post-Newtonian transformation 
of the spatial coordinates contains terms due to the non-perturbative influence of the gravitational 



field. This non-Lotentzian behavior of the post-Newtonian transformations was discussed first by [39| 
for the case of post-Galilean transformations. Our derivations differ from the latter by taking into 
account the acceleration of the proper reference frame and by including the infinitesimal precession of 
the coordinate axes with the angular velocity tensor ^2^^, given by expression (34c) as follows: 

b^M ybMo ybMo 

+(7+i)^^]r(^S'+^f')]- (41) 

where the summation is performed over the other bodies of the solar system. This expression re- 
derives and generalizes the result for the precession of the spin of a gyroscope sq attached to a test 
body orbiting a gravitating primary. Previously this result was obtained from the theory of Fermi- 
Walker transport [|13|. Indeed, in accord with Eq.(346), this spin (or coordinate axes of a proper 
Hermean reference frame) will precess with respect to a distant standard of rest such as quasars or 
distant galaxies. The motion of the spin vector of a gyroscope can be described by the relation: 

^ = [Qm X so]. (42) 
at 

By keeping the leading contributions only and neglecting the influence of the Mercury's intrinsic spin 
moment, we obtain from the expression (41) the angular velocity il.M in the following form: 

= (7 + h-^[RM X VMo] - (7 + l)|f [^M X Ve] + 

+(7 + - 3(5;iVM)iVAf), (43) 

where vmo ^-^id Vq are Mercury's and the Sun's barycentric orbital velocities and 5*0 is the solar 
intrinsic spin moment. 

The first term in Eq.(43) is known as geodetic precession ||8^. This term arises in any non- 
homogeneous gravitational field because of the parallel transport of a direction defined by sq. It can 
be viewed as spin precession caused by a coupling between the particle velocity vmq and the static 
part of the space-time geometry. For Mercury orbiting the Sun this precession has the form: 

^g = {i+\)^{Rmxvm,). (44) 

This effect could be studied for the Mercury Orbiter, which, being placed in orbit around Mercury is 
in effect a gyroscope orbiting the Sun. Thus, if we introduce the angular momentum per unit mass, 
L = Rm X ijmo) of Mercury in solar orbit, the equation (44) shows that VLq is directed along the pole 
of the ecliptic, in the direction of L. The vector VLg has a constant part 

j1o = ^(1 + 27)^ = ^^-0.205 7yr, (45a) 

with a significant correction due to the eccentricity cm of the Mercury's orbit, 

cosoJA/t = -(1 + 27) ^® ^^ ej\/ cosc<jj\/to = — ^ • 0.126 cos c^A/^o "/yi'i (456) 
2 aM 3 
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where ujm is Mercury's siderial frequency, io is reckoned from a perihelion passage; qm is the semimajor 
axis of Mercury's orbit. 

Geodetic precession has been studied for the motion of lunar perigee and its existence was first 
confirmed with an accuracy of 10% [37|. Two other groups have analyzed the lunar laser-ranging data 
more completely to estimate the deviation of the lunar orbit from the predictions of general relativity 
(||88[|, ||89[| ). Geodetic precession has been confirmed within a standard deviation of 2%. The precession 
of the orbital plane proposed for ESA's Mercury arbiter (periherm at 400 km altitude, apherm at 16,800 
km, period 13.45 hr and latitude of periherm at +30 deg) would include a contribution of order 0.205 
"/yr from the geodetic precession. We recommend that this precession be included in future studies 
of the Mercury Orbiter mission. 

The third term in expression (43) is known as Lense-Thirring precession H-lt- This term gives 
the relativistic precession of the gyroscope's spin sq caused by the intrinsic angular momentum S of 
the central body. This effect is responsible for a small perturbation in the orbits of artificial satellites 
around the Earth (|^^, 0])- However, our preliminary studies indicate that this effect is so small for 
the satellite's orbit around Mercury that will be masked by uncertainties in the orbit's inclination. 



6 Discussion 

The use of tracking data from orbiters for relativity tests requires the consideration of more error 
sources than for landers. This is because of the need to convert from the measured earth-spacecraft 
distance to the desired earth-planet distance. This involves determining the orbit of the spacecraft 
about the planetary center of mass, which requires solving from the tracking data for a number of 
spatial harmonics of the gravitational field, solving for radiation pressure, and other non-gravitational 
forces. Firing of attitude control jets which produce unbalanced forces are of particular concern. The 
orbit determination of the Mars orbiter Mariner 9, for example, was substantially affected not only 
by these factors, but also by the fact that the spacecraft was placed on a 12 hr period orbit with low 
periapsis Thus, in order to precisely describe the motion of the Mercury Orbiter relative to 

Earth, one must solve two problems, namely: (i) the problem of the satellite motion about Mercury's 
center of mass in the Hermean-centric frame, and (ii) the relative motion of the both planets - Earth 
and Mercury - in the solar system barycentric reference frame. Our analysis provides a framework for 
the complete solution of these two problems in terms of the corresponding differential equations. 

The formalism presented in this paper addresses the general problem of radio tracking of a Mercury 
Orbiter. We have presented the Hermean-centric equations of the satellite motion, the barycentric 
equations of the planet's motion in the solar system barycentric reference frame, and the coordinate 
transformations which link these equations together. In particular, our analysis shows that in a 
proper Hermean-centric reference frame, the equations of the satellite motion depend on Mercury's 
gravitational field only. This set of equations is available in the form of the motion of test bodies 
in the isolated gravitational one body problem. The existence of the external gravitational field 
manifests itself in the form of the usual tidal forces, but it also determines the dynamic properties of 
the constructed Hermean-centric proper reference frame. Within the accuracy expected for the future 
Mercury Orbiter mission, one can completely neglect the post-Newtonian tidal terms. However, while 
constructing this reference frame, we derived terms smaller than the expected accuracy of future 
experiments. Indeed, the last term in equation (33) is due to the coordinate transformation to the 
Fermi-normal-like reference frame in the planet's vicinity. One can neglect this term for solar system 
motion. However, if one applies the results to problems of motion within a more intensive gravitational 
environment, this term can play a significant role. The application of the results obtained here to 
problems of motion of double pulsars is currently under study and will be reported later. 
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Appendix A: Generalized Gravitational Potentials. 

The generalized gravitational potentials for the non-radiative problems in the WFSMA are given as in 



w-{z^) = J d'z'piz'nv.iz'n ^'^ p}''^^^ 

U"P{zP) = f dh'p{z'P)- 



{z'^ - z"^){z'^ - z'l^) 



m{zP) = -(7 + l)$i - (37 + 1 - 2(5)^2 - $3 - 37^4, 
where the other potentials are given as follows: 

^ ' J \z''-z"'\ ^ ' J ' 

In order to indicate the functional dependence in the potentials introduced above, we have used the 
following notation: (z^) = {z^,z'^). Then for any function / one will have: /(z^) = f{z^^z'^) and 
f{z'r>) = f{z^z'n- 

Appendix B: Simplified Barycentric Equations of Motion. 

In this Appendix we will present the simplified barycentric equations of motion corresponding the 
Lagrangian Eq.(6). Assuming that bodies in the system possess the lowest intrinsic multipole moments 
only, one can obtain the corresponding simplified equations of motion. Thus, with the help of the 
expressions (6), for an arbitrary body (a) these equations will read as follows: 



b^a ah bjta ab 



A a I ^afe I ab 
■^ab+ — + -JT' 
'ab I ab 



(2/3 + 27 - 2r + l)ma + {2p + 27 - 2T)mb 

Tab ^ 
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where, in order to account for the influence of the gravitational binding energy Eb, we have introduced 
the passive gravitational rest mass M-b (H^, as follows 

(52) 



Mb = mb{l + {3 + -f- 4p)Eb + 0{c- 



The unit vector Uab must also be corrected using the gravitational binding energy and the tensor of 
the quadrupole moment /"^ of the body (a) under question: 



JX^l jal3 

Kb = <b(l + (3 + 7 - 4/?)^a + 5nabxnab^,^) + 2n„b^^ + ©(c" 
The term A'^b in the expression (Bl) is the orbital term which is given as follows: 



+Kb{'^a\Va - (7 + 1 + r)VabxVab " ^T{nabxVab^f - ^{nabxV^f 

The spin-orbital term B"f^ has the form: 



ab 



(^ + 2^)vabxiS^' + S^') + IvaxiS:' - S^') + 



m 



+:^{l + 2j)nabxVabp[n'',biS^' + Sb 



'aX\ 



na6"(5f + 5f) 



+ 



+^nabx 



nabiVapS^^ - Vb^S^^) + Uab^Vab'^S'l 



aX 
b 



The term C^b is caused by the oblateness of the bodies in the system: 



Cb = '^n-abplb^ + f>n%nabxnabf,Ii 



Xfi 



iB5) 



(B6) 



And, finally, the contribution I'^^c the equations of motion Eq.(i?l) of body (a), caused by the 
interaction of the other planets (6 / a, c / a, b) with each other is presented as: 



77 

-T-)a ab 

^abc — 2 
'ab 



\l-2P)—-2{P + j)—' 

^bc ^ac 



+ 



+T—f-{nbcX+ncax)+T 2 



riabxKc , 1/. , o n^6caA°^ , , nft, 



+ 2(1 + 7)- 2 > 



(i?7) 



Where A^,^ = r//^^ + n^,n;^, and H^,^ = r?/^- + 3n^^,n:,. 

The metric tensor Eq.(3), the Lagrangian function Eq.(6) and the equations of motion Eqs.(i?l-i?7) 
define the behavior of the celestial bodies in the post-Newtonian approximation. These equations may 
be considerably simplified by taking into account that the leading contribution to these equations is 
the solar gravitational field. Hence they may be used for producing the numerical codes in relativistic 
orbit determination formalisms for planets and satellites ([^, |^^, |91|, |^^) as well as for analyzing 



the gravitational experiments in the solar system ([|l^], |l73[], [|85[). Moreover, taking into account the 
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expected accuracy of the orbit determination for the future Mercury Orbiter mission, one may neglect 
the post-Newtonian perturbations caused by the other planets on the Mercury's orbital motion in the 
terms ^a6''^ab'^ab'^ab' leaving only the solar contribution. The set of the resulting equations now 
defines the properties of the motion of the Mercury in the solar system barycentric reference frame. 



Appendix C: Equations of the Spacecraft Motion. 



In this Appendix we will study the motion of the spacecraft around the planet Mercury. In order 
to obtain the Hermean-centric equations of the satellite motion we will write out the equations of 
geodesies to the required degree of accuracy. For n = a we have: 



ds 



(CI) 



We consider the metric tensor of Riemann space-time to be given by the expressions (30) in this 
case. It is then possible to find the connection components of this space-time needed for subsequent 
computations: 



dU dWa dU^ 



7- 



dy^ dy^ dy^ J 



+ 2(7 + 1) 



dvi 
dy'a 



+ (27 - lK,aa,^ - j6'^ ■ a^ao<^aox - r/"" E [(^i^"),, + 27^a(5i^b)„ 



+ 



dy'a 



= l^t^^ + (7 + l)'^r^^ + 7Wa[>„"l + 



+(7+i)E[(^l.n")„-(5'^^.H,)„' 



y^a+0{\yl\')+0{c-^), 



r.J.^) = 7(^.^ + ^.^-.,.. ^) + 



-2vpJd^d,U,) - 25'^M.Ub) 1 ■ + o{K?) + o{c-'). 

\ I CL > ' Oj- 



(C2a) 



(C26) 



(C2c) 



To reduce the equation of geodesic motion Eq.(d) we shall use these expressions above and the 
definition for the four-vector of velocity in the form: 



dy'a 



|(5oo + 25om^^^ + 5M^V) 



-1/2 



Then by taking into account that d/ds = vPd/dy^ (with the components of the three-dimensional 

^(0) 



velocity vector of the point body denoted as v?L = dy^/dy^) and by using the Newtonian equation of 



motion of a point body as: 



dy'a 



dya 
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we may make the following simplification: 



dU 



n dlnvP „ /dU a ou 



Substituting this relation into the equations of motion (CI), we find the acceleration a°Q^ of the point 
body: 



8TI dV^ 
-(27 + l)-ro) 5^ - 2(7 + 1)^ - 2(7 + 1)^(0), 



Oy'a 



6^a 



<o)4) (2(7 + ^)dxUa + ^(7 + 3) E • Va) 



+ 



+ 



+2(7 + l)i;fo)K^"5/.HA) -(^vH" 



+ 



+: 



(C3) 



By expanding all the potentials in Eq.(C3) in power series of l/ubao ^^'^ retaining terms with 
~ y°'/\ybao\ o^ly to the required accuracy, we then obtain: 

(977 

where the post-Newtonian acceleration (^a^^g) due to the gravitational field of the body (a) only, may 
be given as 

<5aafo) = 2(7 + f3)Uad''Ua - (7 + ^)5"*la + (2/? - ^)5"$2a+ 

+(1 - 7)a"$4a + - (27 + l)^;fo)^/.Kr + 7no)Xo)^"^«- 



-2(7 + l)^^fo)<o)^/^^a - 2(7 + l)v(o)^ 



a a 



2(7 + 1) / d^':pav>'j^;^ - ^(47 + 3) / d^':^^^+ 

Ja \!Ja ~ ya\ ^ Jo. \Ua ~ Ua I 



Pad^'Ua 



{ya-ya){y^a-y'an 



\y''a-y'^\ 
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The term 5afeO(()) acceleration due to the interaction of the gravitational field of the extended 

body (a) with the external gravitation in the N-body system: 



+ 



Vbao y ybao 



+^ ( (2/3 + 7 - 1)5° + 2(3/? + 7 - mba,N,aoM+ 



'10 



+3P— — \ya\nn [25, N, 



ybao 



baox + KaoirieX + 5Nf,ao .N^ao x, 



+ Oi\y:\') + 0{c-''). 



{C5b) 



Note that the combination of the post-Newtonian parameters in the first term of the expression (C56) 
differs from that for the well known Nordtvedt effect (|13], |2^]). This may provide an independent 
test for the parameters involved. The reason that our third term in this expression differs from the 



analogous term derived in |59| is that, in order to obtain this result Eqs.(C4-C5), we used the consistent 
definitions for the conserved mass density in the proper reference frame. Moreover, in constructing the 
Fermi normal coordinates previous authors used the incomplete expressions for the spatial coordinate 
transformations, which are differing from Eq.(31a). Note that if one decide to use our definitions, 
the result cited above will take the form of Eq.(C56). The next term (Jfea^'o) post-Newtonian 
acceleration caused by the other bodies in the system on the orbit of the body (a) (the effect of the 
post-Newtonian tidal forces): 

^ft«fo) = E f - ( ^(147 - 4/3 - 7)5^ + 3iV,-„iV,„„^i(427 - 16/3 - 17)) + 
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+2(7 + i)^fo)<)t bfcao/^r?- - c/^]) ) + o{K?) + o{c 



{Che) 



Finally, the last term in the expression (C4) 5;,cfl(o) contribution to the equation of motion of 

the non-linear gravitational interaction of the external bodies with each other given as follows: 
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+(2/3 - l)Ka,N,ao^ - K,^N,,J + [(2/3 - 37 + hv^+ 



+ 



_^ mbrric 



ybao ycbo 



1-^ A 



K^baoxKbo + l^tcNbaoxNX- 



-i7 + lmLNci.o, + Nbao,N^+{2p-h-^^vA +0{\y:\') + 0{c-% {CM) 

^ ^ ybagycbo ) 

Thus, the equations presented in this Appendix are represent the motion of a test body in the 
Fermi-normal-Uke coordinates chosen in the proper reference frame of a body (a). Together with the 
coordinate transformations Eq.(22) this is the general solution of the gravitational N-body problem. 

Appendix D: Astrophysical Parameters Used in the Paper. 

In this Appendix we present the astrophysical parameters used in the calculations of the gravitational 
effects for the Mercury Orbiter mission in Section V: 

Solar radius : Rq = 695, 980 km. 

Solar gravitational constant : /i© = c~'^GMq = 1.4766 km. 
Solar quadrupole coefficient (Brown et al, 1989) : J20 = (1.7 it 0.17) x 10~^, 
Solar rotation period : Tq = 25.36 days, 
Mercury's mean distance : gm = 0.3870984 AU = 57.91 x 10^ km. 
Mercury's radius : Rm = 2, 439 km, 
Mercury's gravitational constant : /xm = c~^GMm = 1.695 x lO"^//©, 
Mercury's sidereal period : Tm = 0.241 yr = 87.96 days, 
Mercury's rotational period : tm = 59.7 days , 
Eccentricity of Mercury's orbit : cm = 0.20561421, 
Jupiter's gravitational constant : nj = 9.547 x lO~^/Lt0, 
Jupiter's sidereal period : Tj = 11.865 yr. 
Astronomical Unit : AU = 1.49597892(1) x 10^^ cm. 
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